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ABSTRACT 



Progress in understanding the formation and evolution of planetary nebulae (PN) has been restricted by a paucity of well-determined 
central star masses. To address this deficiency we aim to (i) significantly increase the number of known eclipsing binary central stars of 
PN (CSPN), and subsequently (ii) directly obtain their masses and absolute dimensions by combining their light-curve parameters with 
planned radial velocity data. Using photometric data from the third phase of the Optical Gravitational Lensing Experiment (OGLE) 
we have searched for periodic variability in a large sample of PN towards the Galactic Bulge using Fourier and phase-dispersion 
minimisation techniques. Among some dozen periodically variable CSPN found, we report here on three new eclipsing binaries: M 
3-16, H 2-29 and M 2-19. We present images, confirmatory spectroscopy and light-curves of the systems. 
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1. Introduction 

The key parameter for any star is its mass. The least model- 
dependent method to obtain masses is via Keplerian orbits in 
binary systems. Eclipsing binaries provide critical information 
necessary to secure the masses and other fundamental stellar pa- 
rameters. The mass is obtained from radial velocity (RV) orbits 
of both stars combined with the orbital inclination from photo- 
metric eclipses. 

Yet directly derived masses have not been fully exploited in 
the case of central stars of planetary nebulae (CSPN), partly 
because so few eclipsing binary systems are known. Indirect 
methods, such as nebula modelling (e.g. Gesicki et al. 2006) or 
model atmosphere fitting (e.g. Napiwotzki 1999), offer a poor 
alternative to relatively model-independent techniques such as 
Keplerian orbits. 

Obtaining direct masses for a sufficiently large number of 
CSPN is key to understanding the link between the AGB phase 
of stellar evolution and the final white dwarf stage. Such a 
study will also greatly improve our understanding of PN shaping 
mechanisms (e.g. Zijlstra 2007) and can provide compelling can- 
didates for SN la progenitors (Tovmassian et al. 2008). Because 
most stars of low-intermediate mass experience the PN stage, PN 
are also major contributors to ISM enrichment and the chemical 
evolution of galaxies. 

De Marco et al. (2008, and references therein) summarised 
the status of photometrically variable CSPN, including 12 close 
binaries with known periods. Three of these are eclipsing: UU 
Sge (Pollaco & Bell 1993), V477 Lyr (Pollaco & Bell 1994), 



and BE UMa (Ferguson et al. 1999). All three have viable mass 
determinations using both RV and light-curves, with 0.5-0.7 M 
for the hot 60-120 kK primary and 0.15-0.36 M© for the irradi- 
ated 5-7 kK secondary. Their light-curves display strong irradia- 
tion effects, in addition to typical emission lines from the heated 
hemisphere of the secondary. All these systems are believed to 
be post-CE (common-envelope) binaries, where the secondary 
spirals in through the primary's AGB envelope. The remaining 
stars of De Marco et al. 2008 are non-eclipsing and therefore any 
mass estimates for them are subject to considerable uncertainty. 



Either RV or photometric monitoring surveys may be used 
to find CSPN binaries. However, if the ultimate goal (as in 
this work) is to get the masses from binaries, then it is eas- 
ier and more efficient to start with photometric surveys to pre- 
select eclipsing binaries, for which RV orbits combined with 
photometrically-derived inclinations will secure the masses. 
Indeed, RV monitoring surveys have proven to be difficult with 
many candidates showing RV variability without definitive peri- 
odicity (e.g. De Marco et al. 2004). 



This paper introduces a novel use of extant online photo- 
metric data from the OGLE microlensing survey. We describe 
in Section 2 a search for periodic variability in over 300 PN to- 
wards the Galactic Bulge. In Section 3 we present the discovery 
of three new eclipsing binary CSPN. We discuss the nature of 
these new discoveries in Section 4 and conclude in Section 5. 
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2. CSPN Variability towards the Galactic Bulge 

2.1. OGLE-III 

To address the paucity of binary CSPN we have searched 
for periodic photometric variations in data from the Optical 
Gravitational Lensing Experiment (OGLE; Udalski et al. 1992) 
in its most recent third phase (OGLE-III; e.g. Udalski et al. 
2002). Each phase of the OGLE project improves its temporal 
and areal coverage of the Galactic Bulge and Magellanic Clouds 
over previous phases. The areal coverage of OGLE-III is very 
well-matched to the PN distribution towards the Bulge (e.g. the 
most frequently sampled fields cover \€\ < 5 and -5 < b < -2). 

The OGLE-III survey uses the 1 .3-m Warsaw telescope at the 
Las Campanas Observatory, Chile. It is equipped with a CCD 
mosaic camera with eight 2K x 4K CCDs giving a 35 x 35.5 
arcmin 2 field of view with 0.26" pixels. The survey is conducted 
in the I-band reaching a limiting magnitude of I~20, though 
some V-band reference images are also available. The data span 
approximately six years starting from June 2001 with gaps be- 
tween each Bulge season. 

2.2. Search Method 

Surveys for PN towards the Galactic Bulge have recently been 
considerably enhanced by the Macquarie/ A AO/Strasbourg Ha 
PN catalogues (MASH; Parker et al. 2006, Miszalski et al. 
2008). By effectively doubling the PN population towards the 
Bulge, MASH now enables a large population of Galactic PN to 
be monitored for variability in a homogeneous fashion by one 
photometric survey. 

For each of the over 300 PN that fell in the OGLE-III sur- 
vey fields, we examined OGLE-III I-band finder charts to iden- 
tify central stars for which to extract time- series photometry. 
Approximately 30% of these had CSPN identified unambigu- 
ously with the help of significantly lower nebular contamination 
in the I-band. In more difficult cases we selected a number of 
likely candidates close to the geometric centre of each PN. 

Time series photometry for each identified CSPN was sub- 
sequently retrieved from the OGLE-III database generated by a 
data reduction pipeline based on the difference image analysis 
(DIA) method (Wozniak 2000). The stellar-optimised pipeline 
could not extract some light-curves in cases of nebular contami- 
nation. We examined the photometric data using both the Fourier 
analysis package PERIOD04 (Lenz & Breger 2004) and the 
phase dispersion minimisation (PDM) IRAF task (Stellingwerf 
1978). Initially the data were searched with the Fourier tech- 
nique for frequencies in the range 0-50 d" 1 (the latter being the 
typical Nyquist frequency for well- sampled fields) with a step 
rate 2.2 x 10" 4 d" 1 (~l/5 times the inverse sample length of 
~6 years). Periodicity was identified when the Fourier peak am- 
plitude reached at least 3<x above the noise. The final selection 
of the best frequency was made after visual inspection of the 
phased light-curve. We used the PDM technique to check our 
Fourier results for P~0.01-4 days. 



3. Results 

Some dozen periodically variable CSPN were found during our 
examination of OGLE-III data. A photometric binary fraction of 
10-15% (Bond 2000) compares reasonably well with our results, 
but complex selection criteria involving CSPN identification and 
uneven field sampling precludes a meaningful estimate of our 



binary fraction from being given here. A future paper will be 
dedicated to this topic (Miszalski et al., in prep). 

3.1. Three New Eclipsing Binary CSPN 

Of all the periodically variable CSPN found, those of M 3-16 
(PN G359.1-02.3), H 2-29 (PN G357.6-03.3) and M 2-19 (PN 
G000.2-01.9) stand out as having clear eclipsing light-curves. 
These PN have been known for some time (Acker et al. 1992) 
but are all relatively unstudied (though see Section 13.3b . The 
other variables, including some MASH CSPN, tend to exhibit 
sinusoidal-like light-curves and will be discussed elsewhere. In 
Figure [T] we present for each object an Ha image, the OGLE- 
III I-band image and its OGLE-III light-curve. The data are de- 
scribed in the following sections. 

3.2. Photometry 

3.2.1. Ha Images 

High-quality Ha images of each nebula are presented in the first 
column of Figure [T] We obtained an Ha image of M 3-16 us- 
ing the Mosaic II camera on the CTIO 4-m Blanco telescope in 
June 2008. The H 2-29 and M 2-19 images were taken with ESO 
NTT/EMMI during programs 67.D-05 27(A) and 71.D-0448(A) 
(Ruffle et al. 2004). Note that all images include both adjacent 
[Nil] emission lines in the bandpass of the Ha filters. In all in- 
stances the resolution is well-matched to the OGLE-III resolu- 
tion of 0.26'Vpixel. 

3.2.2. CSPN Identification 

Identification of central stars in the OGLE-III I-band images 
requires careful scrutiny and a full discussion of selection ef- 
fects will be given in Miszalski et al. (in prep). Nevertheless, 
in each of our three eclipsing binary CSPN we have unambigu- 
ously identified the CSPN, as indicated by the markers in Figure 
[T] In the case of M 3-16 and M 2-19 we see faint nebulae on 
the I-band images with centrally located stars. Although no I- 
band nebular emission can be seen for H 2-29, the star marked 
is centrally located with respect to the nebula in the Ha image. 
Furthermore, we note that the CSPN marked is brighter in the 
OGLE-III V- than the I-band image, strongly suggesting that the 
star is blue and is therefore the actual CSPN. 

3.2.3. Light Curves 

Phased light-curves of each CSPN are given in Figure [TJ with 
periods and ephemerides given in Table [B No significant change 
in period was detected between each individual Bulge season for 
each object. 

We visually inspected the light-curves phased with the pe- 
riod recovered by the Fourier method. We doubled the period 
of M 3-16 to separate the overlapping, unequal eclipses. The 
slightly rounded shape of its light-curve between minima is in- 
terpreted as ellipsoidal variations. The light-curve of H 2-29 
phased with the initial Fourier period shows a moderately weak 
secondary minimum at 0.5 phase leading us to accept it as the 
most likely period. However, if the secondary minimum turns 
out to be spurious, then the true period would likely be twice 
our value. When phased with the Fourier period, the light-curve 
of M 2-19 lacks a convincing secondary minimum at 0.5 phase 
and shows an unusually wide primary minimum. Therefore, we 
doubled the Fourier period of M 2-19. 
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Fig. 1. A montage of the three newly discovered eclipsing systems. Each row depicts from left to right: An Ha image of the nebula, 
an OGLE-III I-band image, and a phased OGLE-III light-curve of the central star (marked on the I-band image). Periods and 
ephemerides are given in Table Q] All images are 30 x 30 arcsec 2 with North to top and East to left. The solid curve is a binned 
light-curve representation sampled at 5(p = 0.01. 



3.3. Spectroscopy 

Although nebular spectroscopy exists for all our objects, no 
identification of CSPN features have been reported in the litera- 
ture (Cuisinier et al. 2000; Exter et al. 2004). Spectroscopic ob- 
servations of M 3-16 and M 2-19 were made with the AAOmega 
system on the 3.9-m Anglo- Australian Telescope (AAT), on the 
nights of 16 March 2008 and 29 May 2008. AAOmega is an opti- 
cal multi-object spectrograph (Sharp et al. 2006) that is fibre-fed 
by the 2dF fibre positioner (Lewis et al. 2002). The data were 
taken in service mode and exposure times were 2 x 1500s for 
both M 3-16 and M 2-19, with an additional 1800s on M 2-19. 
We used the 580V and 385R volume-phase holographic (VPH) 
gratings for the blue and red arms, respectively. This gave a cen- 
tral resolution of 3.5 A in the blue and 5.3 A in the red and 
a wavelength coverage of 3700-5800 A and 5720-8850 A, re- 
spectively. 

The two epochs for M 3-16 and three epochs for M 2-19 were 
summed to give a S/N in the continuum at 5300 A of -40 and 
-30, respectively. Both spectra exhibit stellar Hell absorption 
lines expected for hot O-type central stars (Figure [2]). Evidence 
for stellar lines that betray a cool secondary component is weak. 
The Call triplet feature is not seen and the possible presence of 
very faint N III and C III emission lines at -4640 A requires 
confirmation (see Pollacco & Bell 1993, 1994). These features 
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Fig. 2. AAOmega spectroscopy for M 3-16 (top) and M 2-19 
(bottom) showing stellar Hell absorption lines (dotted lines). 
Contaminating nebular emission lines are also present. The spec- 
tra have been rectified and smoothed. 



are more likely to be present in H 2-29, but no sufficiently deep 
spectra currently exist. Our confirmatory spectra of low S/N (in- 
dividually) and low resolution were not intended to uncover ra- 
dial velocity variations. Nevertheless, despite large uncertainties, 
we noticed no variation in the radial velocity of the Hell lines 
with respect to nebular lines. Higher resolution data is needed to 
place limits on any possible variation. 
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Table 1. Basic data for the three new eclipsing binaries. An average I-band magnitude fhi is calculated from the OGLE-III time- 
series photometry of N observations over six Galactic Bulge seasons. The epoch of primary minimum Eq is also given. The current 
zeropoint of fhi is roughly calibrated to an uncertainty of -0.1-0.2 mag. 

PNG Name RA (J2000) Dec (J2000) OGLE Field E (HJD-2450000) P (days) m 7 N 

359.1-02.3 M3-16 17 52 46.1 -30 49 35 BLG180.1 3850.9184 + 0.0015 0.573648 + 0.000020 15.93 1039 

357.6-03.3 H2-29 17 53 16.8 -32 40 38 BLG155.1 3619.5208 + 0.0025 0.244110 + 0.000005 18.12 609 

000.2-01.9 M2-19 17 53 45.6 -29 43 46 BLG101.3 4224.7822 + 0.0035 0.670170 + 0.000020 15.87 1273 



4. Discussion 

4.1. Nature of the Systems 

At this stage we can only comment on the components of each of 
our new systems and leave analysis using the Wilson-Devinney 
method (Wilson & Devinney 1971) to future papers. The pe- 
riods of our systems are short and are within the range of the 
periods known for the majority of binary CSPN (De Marco et 
al. 2008). Our light-curves show none of the very strong irradia- 
tion effects seen for UU Sge, V477 Lyr or BE UMa, suggesting 
a very different system composition, though irradiation effects 
may be present in H 2-29, albeit weaker and uncertain. For M 
2-19 and M 3-16 in particular, the eclipses are of similar and 
modest depths suggesting comparable effective temperatures of 
both components. 

Our low-resolution confirmatory spectroscopy show no de- 
tectable difference between radial velocities of Hell absorption 
lines and nebular lines in M 2-19 and M 3-16. This is most likely 
due to our very limited sensitivity to radial velocity variations of 
up to 50-150 km s" 1 that are seen in known eclipsing CSPN 
(Pollacco & Bell 1993, 1994). If such large variations are not 
confirmed by future monitoring, then we may be seeing the hot, 
ionising component of a triple system, with eclipses generated 
by two main- sequence stars of similar effective temperatures. 
This scenario is supported by the lack of strong irradiation ef- 
fects and the similar eclipse depth in the light-curves of M 2-19 
and M 3-16. However, we do not see any convincing stellar fea- 
tures indicative of this scenario in our low S/N spectra. Such a 
case has been suggested for SuWt 2 that has a similar light-curve 
to M 3-16, but its bright A-type double-lined spectroscopic and 
eclipsing binary CSPN has not yet been ruled out as a foreground 
object (Bond et al. 2000). There may be another explanation and 
more radial velocity data at higher resolution are needed to re- 
solve the issue. 

4.2. Morphology 

The morphologies of our three nebulae all show bipolar features. 
This is consistent with the tendency for close binary stars to eject 
non- spherical nebulae through interacting stellar winds and tidal 
effects (e.g. Soker 1997). Indeed, the morphology of H 2-29 has 
been considered as a likely product of CE evolution with a stellar 
companion (Soker 1997). 

The 'butterfly' shape has been proposed to be the preferred 
morphology of nebulae ejected by close binaries because it is the 
most efficient transfer mechanism of angular momentum to the 
ejecta (Zijlstra 2007). M 2-19 exhibits such a classical 'butterfly' 
morphology which is the first firmly established case with a con- 
siderably short period supporting this hypothesis. Other 'butter- 
fly' nebulae with a binary nucleus include NGC 2346 (Mendez 
et al. 1982), but it has a longer period of 16 days. 

M 3-16 exhibits the most complex morphology of the 
three. Surrounding the south-eastern side of the central star 
we see what may be part of an equatorial torus-like structure. 



Perpendicular to this we have two arms often seen in other bipo- 
lar nebulae extended towards the north-west. If we use this to 
define an axis of symmetry, then 25 degrees away from this axis 
we have two opposing broad jets or ansae. 

5. Conclusions 

We have doubled the known population of eclipsing binary 
CSPN after searching OGLE-III photometry in a large PN cat- 
alogue towards the Galactic Bulge. The post-CE CSPN show 
bipolar morphologies consistent with the current hypothesis that 
close binaries lead to non-spherical nebulae. A following paper 
will present more new close binary CSPN, discuss the selection 
effects of the search in detail and provide a new, independent es- 
timate the binary fraction of PN. Planned high-resolution spec- 
troscopy will refine system parameters and derive masses of our 
new eclipsing binaries. 
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